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ABSTRACT 

We investigate the energy-time lag dependence of the source NGC 4593 using XMM-Afewton/EPIC-pn data. We found that the 
time lag dependency is linear in nature with respect to the logarithm of different energy bands. We also investigate the frequency 
dependent time lags and identify that at some frequency range (5 x 10~ 5 Hz - 2 x 10~ 4 Hz) the X-ray emission is highly coherent, 
mildly frequency dependent and very strongly energy dependent. These observations can be explained in the frame work of the 
thermal Comptonization process and they indicate a truncated accretion disk very close to the black hole. We discuss the plausible 
spectral state to explain the phenomenon and conclude that the observed properties bear a close resemblance to the intermediate 
state or the steep power-law state, found in galactic black hole sources. 

Subject headings: accretion, accretion disks - black hole physics - galaxies: individual: NGC 4593 -galaxies: 
Seyferts. 



1. INTRODUCTION 

It has become increasingly clear in recent times that 
the central engines of Active Galactic Nuclei (AGN) are 
magnified forms of G alactic Black Hole Candidates - GBHCs 
dMcHardv et al. 20061) . The scaled up X-ray timing and 
spectral variability characteristics in AGNs (with respect 
to GBHCs) indicate a common energy generation process 
in these sources. State transition in GBHCs is a common 
phenomenon whereas in AGNs, if scaled with mass, the 
state transition may occur in thousand to millions of years 
depending on the mass of the black hole. A detailed 
temporal and spectral study of a sample of AGNs, how- 
ever, may shed light on the spectral states in AGNs. The 
cross spectrum as well as power -s pectral-density (PSD) 
of Ark 564 (lArevalo et al. 20061: IMcHardv et al. 20071) 
indicate that the source is in very high state and dou- 
ble break in PSD and large separation between breaks 
also argue s against a low/hard state like that seen i n 
Cyg X-l (Pap adakis et al. 20021: iDone & Gierlinski 20051) . 
The very high accretion rate inferred in Ark564 (Ro- 
mano et al. 2004) strongly suggests a Very High State 
in this source. Detailed temporal studies indicate that 
the Seyfert galaxies NGC 3227, N GC 4051 and MGC- 
6-30-15 are in hi gh soft state (IMcHardv et al. 2004t 
IMcHardv etal. 20051: lUttlev & McHardv 20051) . " The 
temporal properties of NGC 3783 and NGC 4258 
suggest that these s ources are in low-har d state 
(Mark owitz etal. 20031: iMarkowitz & Uttlev 2005b . An 
investigation of the PSD of NGC 3783 in a wider frequency 
range, however, showed that the PSD is consistent with a soft 
state model (Summons et al. 2007). Recent investigations 
have revealed that type 1 radio quiet AGN and radio loud 
AGN spectral states are analogous to the very high spectral 
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state of GBHCs (ISobolewska et al. 20091) . 

One of the very interesting and important spec- 
tral states in GBHCs is the very high state which is 
often observed during t he transition from low-hard 
state to high soft state dMcClintock & Remillard 2004t 
iRemillard & McClintock 20061: IDone et al. 20071) . This state 
is characterized by the presence of a steep power-law com- 
ponent (with power-law index T > 2.5) which dominates the 
total observed X-ray flux, and hence this state is also referred 
to as the steep power-law (SPL) state. This state was first 
seen in GX 339-4 and l ater observed in many other GBHCs 
dMivamoto et al. 19931 IDone et al. 2007b . Quasi-periodic 
oscillations in the frequenc y range 0.1-30 Hz is also a com - 
mon property of this state dMcClintock & Remillard 20041) . 
The study of hardness-intensity diagram of various black 
hole sources suggests the onset of re lativistic jet during 
this transition state dBelloni et al. 20051) . A theoretical 
modeling of X-ray spectra in the very high state of the XTE 
J 1550-564 suggests that the disk is truncated and the inner 
part of accretion disk is filled with hot and compact central 
corona (IDone & Kubota 20061) . Anti-correlated hard lags 
were also discovered in XTE J 1550-564 which favors a hot 
compact central co rona in the very high state of this source 
(ISriram et al. 20071) . 

It is very important to identify such states in AGNs. Phe- 
nomenologically, one of the defining characteristics of such 
states is the existence of a confined compact corona, which is 
the source of the hard X-ray po wer-law due to inverse Comp - 
ton scattering of soft photons (Sunyaev & Titarchuk 1980). 
One of the ways to understand such a process is to study 
the time lags between two different energy bands. If the 
basic process is Comptonization, the soft photons gain 
energy because of fewer number of scattering with hot 
electrons present in the corona and as the energy difference 
increases between two energy bands, the time lag also 
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increases. This energy gain time scale have been mea - 
sured in a few A GNs, NGC 7469 (iPaoadakis et al. 2001b . 
MCG6-30-15 (Vaug han et al. 2003^7 NGC 4051 



dMcHardv et al. 2004. NGC 3783 dMarkowitz et al. 20051) 
and Mrk 110 dDasgupta & Rao 2006bT Similar type of 
energy depende nt lag was observe d in GBHC, Cyg X- l 
(ICui et al. 1997b ICrarv et al. 19981: iNowaket al. 1999bh . 
This clearly indicates an underlying common radiative 
phenomenon in both GBHCs and AGNs. But this scenario 
of uniform Comptonization is, however, challenged by the 
study of Fourier dependent lags. The Fourier study suggest 
that the low and high frequency variation may be because 
of different X-ray emitting regions in the accretion disk and 
also can be explained by a Compton cloud of radial depen - 
dence of electron density (Kazan as Hua & Titarchuk 1997b . 
Brenneman et al. (2007) have found a hard delay of the order 
of ~230 s by cross correlating soft and hard X-ray bands 
in NGC 4593. Here we try to investigate this observation 
further by making a detailed energy dependent and frequency 
dependent analysis of the delay. 

NGC 4593 is a barred spiral galaxy of Hubble type SBb, 
which hosts Seyfert 1 nucleus situated towards the constella- 
tion Virgo. It is a close by galaxy at a redshift of z=0.009, 
an apparent visual magnitude of 1 1 .67 and harbors a super- 
massive black hole o f mass of the order of 8.1xl0 6 M Q 
(Gebh ardtet al. 20"00l) - which is consistent with the poorly 
constrained reverberation mapping estimate of < 1.4xl0 7 
M Q (Peterson et al. 2004). NGC 45 93 is highly vari- 
able in X-ray, UV, optical and IR bands (ISantos et al. 1995b 
implying towards a variable continuum source close to the 
black hole. Multi-band spectral energy distribution of NGC 
4593 suggests that the respective accretion disk is trun- 
cated and the inn ermost region is dominated by ADAF 
dLu & Wang 2000b . The X-ray continuum is well fitted with 
photoabsorbed power-law above 2 keV and show spectral 
complexity behavior below this energy. The standard rela- 
tivistic iron line is absent as well as reflection component is 
not seen, instead it has two narrow Fe Kg lines at 6.4 keV 
and 6.97 keV in the X -ray continuum dBrenneman et al. 20071; 
iRevnolds et al. 2 004). Because of missing relativistic spectral 
features in X-ray band, the geometry of NGC 4593 is specu- 
lated to be different from traditional geometry of AGNs or 
relativistic features may be buried in the noise as discussed 
by Reynolds et al. (2004). 

In this paper we make a detailed investigations of temporal 
properties of narrow line Seyfert 1 galaxy NGC 4593 using ~ 
70 ks data obtained by XMM-Newton observatory. 

2. DATA REDUCTION AND ANALYSIS 

NGC 4593 was observed on 2002 June 23 by XMM- 
Newton satellite. We have used EPIC pn camera data during 
which it was operating in small window mode with medium 
filter. We have taken ODF (Observation Data Files) to 
obtain the light curves of the source. Source photons were 
extracted from 40 x 40 arcsec 2 region and equal area is taken 
for collecting background photons, away from the source. 
The epatplot task gave no indication of pile-up during the 
observation. We have selected single and double events 
for the analysis (PATTERN<4 and FLAG=0). The light 
curves are binned for 50 s and XRONOS package is used for 
studying the temporal behavior of the source. 



2. 1. Energy Dependent lag 



We have extracted the light curves in five different energy 
bands (0.3-0.5 keV, 0.5-0.8 keV, 0.8-1.5 keV, 1.5-4.0 keV and 
4.0-10 keV). During the observation the source was highly 
variable and showed two dips separated by a steep increase. 
Initially the average count rate was ~ 12s~' followed by a 
steep increase to ~ 25 .s -1 . This is followed by a dip to around 
~ 18s 1 and settled down at ~ 25i _1 (see Fig. [1). To know 
the time lag in various energy bands, we have used crosscor 
program to evaluate the CCF and the respective delay in the 
various energy bands. We have cross correlated the hard X- 
ray bands with the soft X-ray band (0.3-0.5 keV) and the re- 
sults are shown in Fig. [2] A Gaussian model was used to eval- 
uate the ti me lags and errors are e stimated using the criteria 
A;^ 2 =4 (see Dasgu pta & Rao 20 06) (see Table 1). The mean 
energy is given in Table 1 is the weighted mean for all the pho- 
tons in each energy band. It was found that the magnitude of 
the lag is increasing as the energy bands are becoming harder 
with respect to 0.3-0.5 keV band. It clearly suggests that the 
photons take relatively more time to become relatively more 
harder, most probably due to more number of scatterings (Fig. 
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Fig. 1.— XMM-Newton PN 0.3-10 keV light curve of NGC 4593 binned at 
50 sec. 





Fig. 2. — The cross correlation function (CCF) in various energy bands 
with respect to the soft band (0.30 - 0.50 keV) is shown for NGC 4593. The 
energy bands are 0.50 - 0.80 keV, 0.80 - 1 .50 keV, 1 .50 - 4.0 keV, 4.0 - 10.0 
keV, respectively from the top. The first (top) CCF is vertically shifted up by 
0.3 for clarity. The continuous lines are Gaussian fits and the vertical line at 
zero highlights the observed delay. 



2.2. Frequency Dependent Lag 

If Xi(f) and X2(f) is Fourier transform of two light curves 
in energy bands Ei and E2. Then cross spectrum is defined 
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Fig. 3. — The observed time lags are plotted vs energy for NGC 4593. The 
straight line corresponds to thermal Comptonization model fit (see text). 



Frequency (in Hz) 

Fig. 4. — Coherence function calculated using energy bands 0.3-0.5 keV 
and 0.8-1.5 keV. 



TABLE 1 

Details of the selected energy bands and time lags. 



Correlated energy bands 


mean energy (keV) 


Time lag (sec) 


0.30-0.50 vs 0.50-0.80 keV 


0.64 


171.00*54.0 


0.30-0.50 vs 0.80-1.50 keV 


1.11 


223.43±55.0 


0.30-0.50 vs 1.50-4.0 keV 


2.29 


305.00±58.5 


0.30-0.50 vs 4.0-10.0 keV 


5.87 


652.91±91.0 



as C(f) = X*(f)X2(f). Time lag at frequency f is given by 
6t = arg[C(f)]/2jrf. If S i (f) and S2(f) is source power in two 
different energy bands and Ni(f) and Ni(f) is noise spectrum 
in same bands then coherence function can be defined as 



r(/»- l<COT>l 



• N 



<|Si(/)| 2 )<|S 2 (/)| 2 > 



(1) 



where 



N = [\S l (f)\ 2 \N2(f)\ 2 + \S 2 (f)\ 2 \Nm 2 + \N l (f)\ 2 \N 2 (f)\ 2 ]/m 

(2) 

where m is the number of independent measurements 
(Vaughan & Nowak 1997). The error in coherence is cal- 
culated by the me thod described in dNowak et al. 1999at 
iNowak et al. 199 9b). and the error in the lags are calculated 
using the method described in Bendat & Piersol (1986). To es- 
timate error bars on coherence we used high power and high 
measured coherence limit. We re-binned the cross spectral 
density logarithmically and then computed the coherence and 
lag. In Fig. |4] we show coherence between 0.3-0.5 keV and 
0.8-1.5 keV band light curves as function of frequency. It is 
clear from this figure that coherence drops above 5 x 10~ 4 
Hz. Power spectrum also becomes noisy above this frequency 
(see Fig. |5j. Hence we investigate the nature of lag spec- 
trum in the frequency range 5 x 10~ 5 Hz and 5 x 10~ 4 Hz. 
We also compute cross-spectra between two different energy 
bands 0.3-0.5 keV vs 0.8-1.5 a keV and 0.3-0.5 keV vs 4.0 
keV - 10.0 keV. Frequency dependent lag is shown in Fig. [6] 
and it is evident in the figure that lag scales up with energy 
separation. 

3. DISCUSSION AND CONCLUSION 

3.1. Energy dependent lag from Comptonization 

In this paper we investigated the dependence of time lag 
on energy of observed X-ray radiation in NGC 4593. Our 
results suggests that the time lag increases with the mean en- 
ergy of observed photons. A similar dependence of lag on 
energy has been seen in other AGNs as well and the ob- 
served lag is interprete d as Compton scattering time scale 
(Dasgupta & Rao 20"06h . If this lag is due to the process of 
Comptonization, we can get some useful constraints on the 
geometry of the inner accretion disk. We note here that the 
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Fig. 5.— The Power Density Spectrum (PDS) of NGC 4593 in 0.8-1.5 keV 
energy band. 
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Fig. 6. — Frequency dependent lag between 0.3-0.5 keV vs 0.8-1.5 keV and 
0.3-0.5 keV vs 4.0 keV-10.0 keV energy bands. The line shows the broken 
power-law fit (see text). 

Comptonization model applies to the energy dependent lags 
measured at high temporal frequencies. 

Let R be the size of the region in which the Comptonization 
process is undergoing. If A is the mean free path of the photon 
and tj is optical depth of the plasma in this region, then the 
Comptonization time scale will be given by 



to 



omp 



c * max(l,TT) 



(3) 



Let us assume E a to be the initial energy of the photon and 
after N scattering it attains an energy En (where N= 1,2,3,4), 
i.e. En=Eq*A n , where A=1+4Q+16© 2 (know as amplifica- 
tion factor and &=kT e /m e c 2 ). Then the Comptonization time 
scale can also be written in terms of time lag and N as 



N 



iComp 

where T\ ag is the respective time lag. 



(4) 



The size of the region in which the Comptonization pro- 
cess is undergoing is equivalent to the distance traveled by 
the photons (i.e A, valid only if the medium is optically thin, 
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i.e. max(l, Tt)~1). Using the above two equations we get a 
Comptonizing region size of the order of ~6x 10 12 cm (using 
Ti ag =600 s, for the 0.3 - 0.5 keV and 4-10 keV bands). 

The observed time lags are nothing but product of Comp- 
tonization time scale and difference of successive number of 
scattering, which can be given by 



Tlag 



R 



c*max(l,TT) InA 



(5) 



The observed time lag increases linearly with logarithmic of 
energy (see Fig. |3J. Hence the above equation suggests that 
the Comptonization process is responsible for the observed 
lag. We have fitted the above equation to the observed energy 
dependent time lag (see Fig. [3} and found the size of hard X- 
ray emitting region is ~10 13 cm and when the mass is taken 
7x 10 6 M black hole (Nelson et al. 2004), R comes out to be 

3.2. Disk Geometry and Time Lag spectrum 

The geometry of the accretion dis k in NGC 4593 is quite 
peculiar on the bas is of Fe line dBrenneman et al. 20071: 
iRevnolds et al. 2004b . Also iron line does not show signif- 
icant variation on observed time span and are very narrow 
which are suggestive of absence of optically thick cold matter 
in the central 1000 R g region. The inner 1000 R g region may 
contain radiatively inefficient and hot optically thin plasma. 
The observed time lag suggests that there is also a compact 
central corona with size of ~ 5 R g . The derived Compton 
cloud size from time lag suggests that most of the hard X- 
ray photons are emitted within ~3R g region. It has also been 
argued that NGC 4593 has a truncated accretion disk with a 
two component flow of a geometrically thin disk and an ad- 
vection dominated flow within the truncation radius, a model 
commonly invok ed to explain the l ow-hard and intermediate 
states of GBHs dLu & Wang 2 000). The intermediate state 
state is characterized by a steep power law component and a 
truncated accretion disk (with the truncation radius less than 
that found in the low hard state). Though the power-law in- 
dex is not steeper than that found in other Seyfert galaxies 
dBrenneman et al. 20071) . the available data strongly suggest a 
thin disk converting into an advection dominated central part, 
a configuration similar to the low-hard and intermediate states 
of GBHs. 

Cross spectral analysis reveals that lag observed in this 
source depends upon Fourier frequency. The lag spectrum 
computed using energy bands 0.3-0.5 keV and 0.8-1.5 keV 
reveals a broad peak in the frequency range 5 x 10~ 5 Hz 
to 5 x 10~ 4 Hz (see Fig. |6). In this frequency range, the 
X-ray emission at different energies are highly coherent (see 
Fig. HJ. The lag spectrum, however, shows quite strong 
energy dependency. The structure in the lag spectrum can 
be fitted with a broken power law with first slope a = 0.3 
and second slope ft = -5 with a x 2 = 1-5 for 3 degree of 
freedom (d.o.f), an improvement in^ 2 of ~10 from a fit to a 
constant value. The break is observed at frequency fi, = 2.2 X 
10~ 4 Hz. The observed lag at different frequencies increases 
with energy separation. The lag spectrum of Ark 564 shows 



step-like features (McHardy et al. 2007) and at higher 
frequencies it can a lso be fitted with a broken power-law 
dArevalo et al. 2006!) . In Ark 564 also the lag at different 
frequencies scales up with energy separation. The lags below 
the break frequency is a constant in NGC 4593, which is 
quite different from the trend of lag decresing with increasing 
Fourier frequency seen in other sources like MCG-6-30-25 
(Vaughan, Fabian & Nandra 2003), Mkn766 (Markowitz 
et al. 2007), NGC 4051 (McHardy et al. 2004). The 
fractional lag in this frequency range is 1 - 6%, which is sim- 
ilar to the fractional lag seen in the very high state of Cyg X- 1 . 

If the energy dependent lag is mainly due to Comptoniza- 
tion, then the bulk of the variability (confined to below 2 
x 10~ 4 Hz, see Fig. [6) is nearly frequency independent. It 
is quite possible that there is another variability component 
operating particularly above the break frequency fi, = 2.2 
x 10~ 4 Hz, possibly from accretion disk fluctuations prop- 
agating from outside (see below). The spectral analysis of 
NGC 4593 suggests that the soft excess too can be modelled 
as due to Comptonization from a cooler plasma (Brenneman 
et al. 2007). It is quite possible that there is a hot compact 
central corona within a few Swartzchild radius (responsible 
for the bulk of the radiation and giving rise to the frrequency 
independent delays), surrounded by an optically thin plasma 
cloud extending to several hundred Swartzchild radius 
(responsible for the soft excess). There is a possibility that 
electron density in optically thin corona is varying with 
radius. It was suggested that if corona has non-uniform 
electron densit y then observed frequency dep endence lag can 
be explained dKazanas Hua & Titarchuk 1997b . There is also 
a possibility tha t there exists several dif ferent patches of hot 
Compton cloud (Sob olewska et al. 2 004) giving rise different 
time lags at different time scales. 

Another promising mod el which explains t he observed lag- 
spectrum in Mrk 33 5 dArevalo et al. 20081) is fluctuation s 
propagation model dLvubarskii 19971: iKotov et al. 200 lb . 
which can also give energy dependent delays. This model as- 
sumes that fluctuations in the accretion flow are produced at 
different radial distances from putative black hole. These fluc- 
tuations propagates on viscous time scales and also assumes 
that spectra hardens toward the center. Kotov et al. (2001) has 
computed phase lags for different scenarios which can explain 
the time lags caused by accretion. The lag spectrum observed 
in NGC 4593 is similar to that obtained for freely falling blob 
model(see Fig. 9 of Kotov et al. 2001), but the strong energy 
dependency, however, favors a Comptonization model. 
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